Kondo Breakdown and Hybridization Fluctuations in the Kondo-Heisenberg Lattice 
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We study the deconfined quantum critical point of the Kondo-Heisenberg lattice in three dimen- 
sions using a fermionic representation for the localized spins. The mean-field phase diagram exhibits 
a zero temperature quantum critical point separating a spin liquid phase where the hybridization 
vanishes and a Kondo phase where it does not. Two solutions can be stabilized in the Kondo phase, 
namely a uniform hybridization when the band masses of the conduction electrons and the spinons 
have the same sign, and a modulated one when they have opposite sign. For the uniform case, we 
show that above a very small temperature scale, the critical fluctuations associated with the van- 
ishing hybridization have dynamical exponent z = 3, giving rise to a resistivity that has a TlogT 
behavior. We also find that the specific heat coefficient diverges logarithmically in temperature, as 
observed in a number of heavy fermion metals. 

PACS numbers: 71.27.+a, 72.15.Qm, 75.20.Hr, 75. 30. Mb 



A large number of experiments have been performed 
on metallic heavy fermion compounds close to a zero 
temperature phase transition (a quantum critical point 
(QCP)) driven by applied magnetic field, chemical dop- 
ing or pressure I n the quantum critical regime, the 
thermodynamics and transport properties are very un- 
usual, violating the predictions of the Landau Fermi liq- 
uid theory of metals. The resistivity is quasi-linear in 
temperature over several decades, and in many cases the 
specific heat coefficient diverges logarithmically as the 
temperature is decreased. These unusual observations 
have motivated many theoretical studies that have at- 
tempted to capture these effects. Most theories [J El 01 
are based on the assumption that at the QCP, the Fermi 
liquid is destabilized by spin density wave formation, and 
therefore the critical fluctuations are magnetic in nature. 
In d — 3, these theories fail to capture simultaneously the 
linear temperature dependence of the resistivity and the 
divergence of the specific heat coefficient at low tempera- 
tures |fj . More recently the problem has been approached 
from another perspective which takes the point of view 
that at the QCP, magnetic fluctuations suppress the for- 
mation of the heavy Fermi liquid, driving the effective 
Kondo temperature of the lattice (Tk) to zero 0, 0> IE @ • 
In this picture, the QCP is a bi-critical point where the 
metal experiences fluctuations due to the vanishing en- 
ergy scale Tk as well as the paramagnons. One feature 
that distinguishes between these two classes of theories 
is that the first predicts the Fermi volume to change 
smoothly across the QCP, while the second predicts an 
abrupt change 0- 

In this Letter we explore the possibility that in the 
quantum critical regime, the unusual behavior in ther- 
modynamics and transport is due to critical fluctuations, 
but of a non-magnetic order parameter associated with 
the vanishing energy scale Tk, and not due to param- 
agnons. The order parameter we advocate is the field a 
associated with the hybridization between the localized 



spins and the conduction electrons 

[H 03. At the QCP, 
the effective Kondo temperature for the lattice goes to 
zero, leading to a "Kondo breakdown" of the heavy Fermi 
liquid. The critical fluctuations of a are gapless excita- 
tions, and we study how these fluctuations influence the 
properties of the metal using the formalism of the large 
N Kondo-Heisenberg model. 

Beyond the mean-field level, the Kondo-Heisenberg 
model can be treated as a lattice gauge theory. Senthil 
et. al. ;8| have examined the effect of the gauge fluctua- 
tions in this model, while Coleman et. al. |l| studied the 
zero temperature transport anomalies. In our work, we 
find a number of novel effects associated with the fluctu- 
ations of the a field which were not discovered in these 
earlier studies. 

At the Kondo breakdown QCP where (a) = 0, we ob- 
serve two new phenomena: (1) a can order at a finite 
wavevector leading to spatial modulations of the Kondo 
hybridization analogous to the LOFF state of supercon- 
ductivity (2) the presence of multiple energy scales, 
spread over a very large range in energy. The lowest scale 
is extremely small (of order 1 mK), above which, up to an 
ultraviolet cutoff of order the single ion Kondo temper- 
ature, the critical fluctuations of a exhibit a dynamical 
exponent z = 3. This gives rise to a marginal Fermi liq- 
uid like behavior for the conduction electrons in d = 3, 
with a resistivity that goes as TlogT. A logarithmic 
dependence is also found for the specific heat coefficient 
from both the gauge and a fluctuations. 

Our starting point is the large TV formulation of the 
Kondo-Heisenberg model, where N denotes the enlarged 
spin symmetry group SU(N). It describes a broad band 
of conduction electrons interacting with localized spins 
through anti-ferromagnetic Kondo coupling Jk > 0. The 
localized spins interact with each other via nearest neigh- 
bor exchange Jh > 0. We work with a representation of 
the localized spins in terms of Abrikosov pseudo fermions 
& = HapH,a°otpU,p, where (a, (3) = (1,N), with the 
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constraint of n/ = N/2 spinons per site i. The interac- 
tions which are quartic fermionic terms can be decoupled 
using Hubbard-Stratonovich fields cpij — > ^ Q fj a fj a f° r 
the Heisenberg exchange, and crj — > fia°ia for the 
Kondo interaction. Following Ref. we assume that in 
d = 3, ip condenses in a uniform spin liquid phase that 
gives dispersion to the spinons, which is an essential in- 
gredient for Kondo breakdown to occur (physically we 
interpret the uniform spin liquid as a mean field descrip- 
tion of the short range magnetic correlations that persist 
when a magnetic ground state is destroyed by quantum 
fluctuations). This gives the Lagrangian 
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(V, the volume of the system, is set to 1). The above La- 
grangian has a local U(l) gauge invariance 0. The La- 
grange multiplier A* (scalar potential) enforces the con- 
straint rif = N/2 per site. Given a state which satisfies 
the above constraint, a single hop of a spinon will vio- 
late it. Consequently only simultaneous opposite hops 
of spinons between two neighboring sites are physically 
allowed. This implies that the local spinon current op- 
erator Jfi = at each site. The gauge fields ay (vector 
potential) associated with the phase of tpy ensure that 
this condition is satisfied. 

There are two important parameters in Eq. ^ First, 
a = tpo/D, which is the ratio of the spinon to the con- 
duction electron bandwidth D (note from Eq. 2] that 
for (7=0, 4>q = Jh)- Second, while the spinon band is 
half-filled due to the constraint (henceforth we assume 
N = 2), the conduction band filling is generic. With- 
out any loss of generality we take the conduction band 
to be less than half filled. This implies that the Fermi 
wavevector of the spinon band, kpo, is different from that 
of the conduction band, kp. We denote the mismatch by 
q* = kpo — hp. In the following we take a and (q* /fop) 
to be small. We identify aD with the single ion Kondo 
scale (T K = De" 1 ^ PoJk ) which is typically 10K in heavy 
fermions. Assuming D ~ 10 4 K, we get a ~ 10 -3 . 

At the mean-field level, the parameters tpo, (A,) and 
(crj) are determined by minimizing the free energy, F. 
The mean-field phase transition between the spin liquid 
state, {<Ti) = 0, and the heavy Fermi liquid state with a 
lattice Kondo temperature Tk ~ , Kpo{ai) 2 , occurs when 
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where n/ c (g,0) is the static electron-spinon (fc) po- 
larization. We solve this equation for two different 
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FIG. 1: (a) q dependence of II/ C (a=0.01 and q* /kir=0.1). 
Note differing scales for the e-e and e-h cases, (b) Quantum 
critical point as a function of Jk and Jh for the e-e (q=0) 
and e-h (q=1.2q*) cases. 



situations, the result of which is depicted in Fig. 1. 
(i) e-e case, where both the bands are taken to be 
electron-like. Linearizing the fermionic dispersions we 
have efe = Vp{k — kp) for the conduction electrons, and 
(° k = avp(k — kp — q*) for the spinons (where k = |k|). 
For linearized dispersions, II/ c (g, 0) turns out to be q- 
independent. Inclusion of the curvature stabilizes a sec- 
ond order phase transition around q — 0, the polarization 



taking the form n /c («,0) = ^(Ina + ^t 2 "^" ^) 
where po = l/D is the conduction electron density of 
states at the Fermi energy. In this case the T = 
phase transition occurs at a critical Kondo coupling of 
J K = l/(po ln(l/ poJh)) 0- (h) e-h case, where the 
conduction band is taken to be electron-like as before, 
while the spinon band is hole-like with a linearized dis- 
persion e° = —avp(k — kp — q*). In this case we find 

n/cM) = ^(in -glf+V - 2 + V ln iS^r). which 

has a minimum at q = l.2q* independent of a. In this 
state Tk is modulated, with nodes in space where Tk 
vanishes. This solution is similar to the spin density 
wave instability encountered in chromium |l5j and in the 
LOFF state of superconductivity 0,0]. Fig. la shows 
that for parabolic bands the minimum of the effective po- 
tential is lower in the e-h case than in the e-e case. Thus, 
for parabolic bands, the modulated solution is more sta- 
ble (Fig. lb). However the question of which solution is 
realized in real compounds will be material dependent. 

We now turn to the fluctuations around the mean- 
field solution. We present our results for the simpler e-e 
case, leaving the more complex e-h case for a later paper. 
In the quantum critical regime there are two important 
types of gapless fluctuations, namely the gauge fluctua- 
tions associated with and the critical fluctuations of 
the Kondo bosons a. 

The gauge fluctuations of this theory have been stud- 
ied earlier by Senthil et. at Here we summarize the 
salient points to put our work in perspective. It is conve- 
nient to work in the Coulomb gauge V • a = 0, where the 
vector gauge fields (fx = x,y, z) are purely transverse 
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FIG. 2: (a) Structure of D a for the e-e case (a << 1) at the 
QCP for positive frequencies, z is the dynamical exponent 
in the various regions, which are delineated by the dashed 
lines, equal to avF{q* ±<?) and vp(q" ±<?)- u co u is the zero of 
D~ x in the z=2 regime (a propagating mode), and Y q is the 
dispersion of the damped mode in the z=3 regime (maximum 
of the imaginary part of D a ). Note the presence of energy 
scales E* ~ W~ 4 aD and Ex = chVf q* ~ W^aD. (b) Phase 
diagram for the e-e case (a << 1). Note the crossover from 
z=2 to z—3 behavior at E* (dotted line). The inset shows the 
phase diagram on a linear scale. The solid line is the crossover 
line in the Kondo phase, the dashed line the crossover line in 
the spin liquid phase. 



|l3|. The fluctuations of A decouple from those of <z M , 
and give rise to a screened Coulomb interaction. As such, 
they are massive and can be neglected. The gauge fields 
a^, which act as vectorial Lagrange multipliers to ensure 
that the local spinon current Jfi = 0, do not have any 
intrinsic dynamics of their own. Their dynamics is gen- 
erated via coupling with the spinon band, and therefore 
they are over-damped. For frequencies smaller than the 
spinon bandwidth aD, the transverse gauge field prop- 
agator D MI/ (x, t) = (T T [a M (x, -r)a„(0, 0)]) has the stan- 
dard form D f _ l „(q,in n ) = (5^ - q fl q^/q 2 )Il^ 1 (q,iQ, n ), 
with H(q, ifl n ) °c [(q/2kp n ) 2 + \Q n \/(avpq)]. This is the 
typical form for excitations with dynamical exponent z = 
3, which in d = 3 are known jl7j to give a contribution to 
the specific heat coefficient 7 = —d 2 F/dT 2 oc \n(aD/T) 
and to the static spin susceptibility Sxs oc T 2 \n(aD /T) . 
Finally, when the compact nature of the U(l) gauge 
group on the lattice is taken into account, the gauge fluc- 
tuations convert the finite temperature mean-field tran- 
sition line into a crossover line . 

In the quantum critical regime, the fluctuations of 
the complex order parameter fields (<7?,<7i) are mass- 
less as well (ignoring for the moment the Independent 
mass generated by the quartic |ct| 4 coupling). The prop- 
agator for these fluctuations is defined as D cr (x,r) — 
(T T [a^(x,T)a(0,0)]) with D-\q,m n ) = l/J K + 
Uf c (q, 0) + AIT /c (q, ifi„) where 

att / x \- Po[TX 1± \n(X 1± ) ± X 2± ln(X 2± )} 
An /C (g, zO„) = ^ 2avMi a) 

(3) 

with Xi± = —aiQ n ± avpq — avFq*, X 2 ± = —i£l n ± 



avpq — avFq* , and II/ c (0,0) = po ln(a)/(l — a) (AII/ C is 
the dynamical part of the fc polarization). The different 
regimes of D a {q 1 iQ n ) with their associated dynamical ex- 
ponents z are summarized in Fig. 2a. At high energies, 
one finds z = 00 behavior consistent with quasi-local be- 
havior. But we find physical properties are dominated 
by the z — 3 and 2 = 2 regimes. These two regimes can 
be understood as follows. Due to the mismatch between 
the two Fermi surfaces, a minimum momentum is nec- 
essary to excite inter- band (fc) particle- hole pairs. As a 
result for CI < avp{q* — q), excitations of a do not decay 
into particle-hole pairs but propagate ballistically with 
D- 1 {q,m n ) « p [q 2 /(4k 2 F )-in n /E x l i.e., z = 2 (where 
Ex = avpq*)- This behavior is cutoff for frequencies 
tt > E* with 



E* = caD(q* /k F ) 3 and 
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above which the dynamical exponent z changes from 2 
to 3. For most of the phase space, the spectrum for 
the fluctuations of a lie within the inter-band particle- 
hole continuum, making their dynamics over-damped 
with D- 1 (q,ifl n ) sa p [q 2 /{4k 2 F ) + n\n n \/(2av F q)], i.e., 
z = 3. The energy scale in the z=3 regime is T q /Ex ~ 
q 3 / (2irk F q*) which has an infrared cutoff at E* because 
of the mismatch vector q*. The ultraviolet cutoff scale 
for the z — 3 regime is avp (q + q*) which is of order aD 
for q ~ kp. The energies E* (infrared) and aD (ultravi- 
olet) appear as crossover scales for any physical property 
that is affected by the excitations of a. For a one impu- 
rity Kondo scale aD ~ 10K, and q* jkp ~ 10 _1 0, we 
estimate E* ~ ImK. E* is therefore a very small energy 
scale which is essentially unobservable. 

The crossover lines in T that define the quantum crit- 
ical region are symmetric around the QCP 5 = 5 C , where 
S = l/(p J K ) (Fig. 2b). These are determined by the In- 
dependent mass generated by the quartic |cr| 4 coupling. 
For T < E* , we find that the leading contribution is 
from the z = 2 regime (proportional to T 3 / 2 for d = 3), 
so that the crossover temperature T oc \S — S c \ 2 ^ 3 , while 
for T > E* the z — 3 regime dominates giving a crossover 
temperature T oc \S — S c \ 3 ^ 4 for d = 3. 

Next we examine the contribution to the free energy 
from the fluctuations of a. We find that the leading con- 
tribution is entirely due to the z = 3 regime since it comes 
from a much larger phase space volume (the z = 2 con- 
tribution is similar to that of a gapless magnon mode). 
For E* < T < aD, we find F(T) - / n B Im ln(Z)- 1 ) ~ 
-k 3 F T 2 / (%aD)\n(aD /T), which is a typical result for 
z = 3 excitations. This implies a contribution to the 
specific heat coefficient 7 ~ 2k 3 F / (9aD) \n(aD/T), which 
adds to a similar contribution from the transverse gauge 
fluctuations. For T < E* , the infrared cutoff sets in, and 
the specific heat coefficient from the a fluctuations satu- 
rates. This regime is then dominated by the logarithmic 
contribution from the transverse gauge fields 17]. 
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We now calculate the self energy of the conduction elec- 
trons due to the hybridization fluctuations (c ^ / + a). 
This is defined as £ c (fc, iui n ) = TJ^^ q G 'f(k + q, iuj n + 

ifl n )D a (q,ifl n ), where G7 1 (k,iu> n ) = {iu n — e°) is the 
inverse propagator of the spinons. As in the case of the 
free energy, we find that the leading contribution is due 
to the z = 3 regime of D a (q,iD, n ). At T — and for 
E* < uj < aD, we find ImY, c (kp ,uj) ~ k 2 F j (6iravpp )w. 
The temperature dependence of IxaE c (kp, to — 0,T) 
involves a frequency integral weighted by the factor 
iib + rip = 1/ sinh(f2/T). This makes the integral loga- 
rithmically divergent in the infrared, which is cutoff by 
E*. For E* <T < aD we find 

ImE c (fc F ,w = 0,T) - k%/(Qnoev F po)T]n(2T/E*). (5) 

For lo,T < E* , the self energy is Fermi liquid like. 

We turn to the T-dependence of the static spin sus- 
ceptibility, Xs(T). At the mean field level, we find what 
is usual for band fermions, namely a constant part plus 
a T 2 term. To calculate the correction beyond mean 
field (Sxs) due to the Kondo bosons, we note that in a 
magnetic field B, there is an additional (B / (aD)) 2 con- 
tribution to D~ x {q,m n ). This gives Sx s (T) oc T 4 / 3 for 
E* <T < aD, and a T 2 dependence below E* (so below 
E* the T 2 ln(T) contribution due to the gauge bosons 
dominates) . 

Finally we discuss the temperature dependence of the 
resistivity, p, that is obtained in the quantum critical 
regime. Eq. gives the T-dependence of the inverse life- 
time t^ 1 oc ImS c (T) of the conduction electrons. For 
one band models experiencing q ~ scattering, this life- 
time cannot be associated with the transport lifetime, 
T tr , because the leading contribution to the self energy 
comes from forward scattering processes which are not 
effective in relaxing the current. Consequently, when 
vertex corrections are taken into account, r^ 1 acquires 
an additional temperature dependence proportional to 
q 2 ~ T 2 / z . However, our model consists of two bands, 
one of light particles (the conduction electrons) which 
scatter from very heavy particles (the spinons) [2fJ. As 
such, the charge neutral spinons act as an effective bath 
for the relaxation of the conduction electron current (the 
other charge carrying modes, the complex a bosons, 
have over-damped dynamics, and therefore the current 
is mostly carried by the conduction electrons). The 
first non-zero vertex correction involves two a boson ex- 
change processes. Such a correction is small by a factor 
of a. Therefore r tr can be identified with t c , and for 
E* < T < aD we find 

Sp(T) = p(T) -p(0) oeTln(2T/£*). (6) 

For T < E*, 5p(T) oc T 2 , but E* is extremely small 
(~ 1 mK). Thus, the Kondo-Heisenberg model captures 
one of the most mysterious features of quantum critical- 
ity in heavy fermion compounds, namely the quasi-linear 



resistivity observed for most compounds over a large tem- 
perature range. 

In conclusion, we studied the Kondo breakdown QCP 
of the Kondo-Heisenberg model in d = 3. Over a large 
temperature range, we find that the critical fluctua- 
tions have a dynamical exponent z = 3, giving rise to 
marginal Fermi liquid behavior for the conduction elec- 
trons. The specific heat coefficient has a log(l/T) de- 
pendence, while the resistivity has a TlogT behavior. 
The Kondo-Heisenberg model is characterized by mul- 
tiple energy scales, and as such shows great promise in 
explaining the various subtleties associated with heavy 
fermion quantum critical behavior. 
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